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(54) Methods and apparatus for reducing z-axis non-uniformity artifacts 



(57) Methods and apparatus of correcting image da- 
ta for z-axis non-uniformity in a computed tomography 
(CT) system (10) are described. In one embodiment, the 
CT system (1 0) includes a configurable multislice detec- 
tor array (1 8) and a calibration algorithm. In accordance 
with one embodiment of the algorithm, output signals 



from a plurality of detector cells are sampled to generate 
an x-ray beam z-axis non-uniformity profile. The x-ray 
beam z-axis profile is then used to determine a detector 
z-axis non-uniformity profile for the detector array (18). 
Utilizing the z-axis profiles, the output signals from the 
detector cells may be calibrated to reduce z-axis non- 
uniformity errors. 
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Description 

[0001] This invention relates generally to computed 
tomography (CT) imaging and more particularly, to cor- 
rection of image data for z-axis non-uniformity error in- 
troduced into such data by x-ray beam and detector cell 
non-uniformity. 

[0002] In at least one known CT system configuration, 
an x-ray source projects a fan-shaped beam which is 
collimated to lie within an X-Y plane of a Cartesian co- 
ordinate system and generally referred to as the "imag- 
ing plane". The x-ray beam passes through the object 
being imaged, such as a patient. The beam, after being 
attenuated by the object, impinges upon an array of ra- 
diation detectors. The intensity of the attenuated beam 
radiation received at the detector array is dependent up- 
on the attenuation of the x-ray beam by the object. Each 
detector element of the array produces a separate elec- 
trical signal that is a measurement of the beam attenu- 
ation at the detector location. The attenuation measure- 
ments from all the detectors are acquired separately to 
produce a transmission profile. 

[0003] In known third generation CT systems, the x- 
ray source and the detector array are rotated with a gan- 
try within the imaging plane and around the object to be 
imaged so that the angle at which the x-ray beam inter- 
sects the object constantly changes. A group of x-ray 
attenuation measurements, i.e., projection data, from 
the detector array at one gantry angle is referred to as 
a "view". A "scan" of the object comprises a set of views 
made at different gantry angles, or view angles, during 
one revolution of the x-ray source and detector. In an 
axial scan, the projection data is processed to construct 
an image that corresponds to a two dimensional slice 
taken through the object. One method for reconstructing 
an image from a set of projection data is referred to in 
the art as the filtered back projection technique. This 
process converts the attenuation measurements from a 
scan into integers called "CT numbers" or "Hounsfield 
units", which are used to control the brightness of a cor- 
responding pixel on a cathode ray tube display. 
[0004] Multislice CT systems are used to obtain data 
for an increased number of slices during a scan. Known 
multislice systems typically include detectors generally 
known as 3-D detectors. With such 3-D detectors, a plu- 
rality of detector elements form separate channels ar- 
ranged in columns and rows. Each row of detectors 
forms a separate slice. For example, a two-slice detec- 
tor has two rows of detector elements, and a four-slice 
detector has four rows of detector elements. During a 
multislice scan, multiple rows of detector cells are simul- 
taneously impinged by the x-ray beam, and therefore 
data for several slices is obtained. 
[0005] I n a CT system having such a 3-D detector, the 
intensity of detector measurements are derived by com- 
bining, along the z direction, multiple detector outputs. 
These outputs are supplied as inputs to a data acquisi- 
tion system. If the detector outputs to be combined are 



obtained from detectors having different individual 
gains, the combined signal represents a weighted sum 
of the incoming detector signals where the different de- 
tector gains cause different weighting. Know gain cali- 
5 bration algorithms may be used to correct for this differ- 
ent weighting. However, these known calibration algo- 
rithms are unable to correct the signals when the error 
is caused by z-axis non-uniformity of the x-ray beam and 
the detector. As a result, a small variation in the non- 
uniformity of the z-response results in a substantial 
amount of artifacts in the reconstructed image. 
[0006] It would be desirable to correct for z-axis non- 
uniformity of the x-ray beam to reduce image artifacts. 
It also would be desirable to more accurately generate 
an image regardless of the z-axis non-uniformity of the 
detector arrays. It further would be desirable to provide 
such imaging without significantly increasing the cost of 
the system. 

[0007] These and other objects may be attained in a 
system which, in one embodiment, corrects the error in 
projection data resulting from z-axis non-uniformity of 
an x-ray beam and detector array. More particularly and 
in accordance with one embodiment, the present algo- 
rithm determines a z-axis non-uniformity profile for an 
incoming x-ray beam flux by fully sampling several 
channels of the detector array. Based on the fully sam- 
pled channels, a z-axis profile is generated for each of 
the sampled channels. Utilizing the z-axis profiles from 
the fully sampled channels, a z-axis profile is generated 
for the non-fully sampled channels. Specifically, as a re- 
sult of the x-ray beam profile being a slowly varying func- 
tion in the z-direction with respect to the channel 
number the z-axis profile of the fully sampled channels 
may be used to generate z-axis profiles for the adjacent 
non-fully sampled channels. 

[0008] In addition, the present algorithm determines 
a z-axis non-uniformity profile of the detector array by 
estimating the error caused by the gain variations of the 
individual detector array cells. Specifically, utilizing the 
z-axis profiles for the fully sampled channels, a detector 
cell z-axis non-uniformity profile may be generated for 
the detector array so that the output signals of the de- 
tector array may be calibrated to reduce z-axis non-uni- 
formity errors. 

[0009] The above described system enables correc- 
tion of error due to z-axis non-uniformity of the x-ray 
beam and the detector array. Furthermore, the system 
enables correction of varying detector cell gain profiles. 
As a result, artifacts are reduced in a displayed image. 
The system also does not significantly increase the cost 
of the imaging system. 

[0010] Embodiments of the invention will now be de- 
scribed, by way of example, with reference to the ac- 
companying drawings, in which: 

[0011] Figure 1 is a pictorial view of a CT imaging sys- 
tem. 

[0012] Figure 2 is a block schematic diagram of the 
system illustrated in Figure 1. 
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[0013] Figure 3 is a perspective view of a CT system 
detector array. 

[0014] Figure 4 illustrates a fully sampled channel and 
a non-f ully sampled channel of the detector array of Fig- 
ure 3. 

[0015] Figure 5 illustrates detector cell data combin- 
ing in a four-slice mode having 5mm thickness image 
slices. 

[0016] Figure 6 illustrates a z-axis non-uniformity 
curve in accordance with one form of the present inven- 
tion. 

[0017] Figure 7 illustrates two channels of the detec- 
tor array shown in Figure 4. 

[0018] Referring to Figures 1 and 2, a computed tom- 
ograph (CT) imaging system 10 is shown as including a 
gantry 12 representative of a "third generation" CTscan- 
ner. Gantry 12 has an x-ray source 14 that projects a 
beam of x-rays 1 6 toward a detector array 1 8 on the op- 
posite side of gantry 12. Detector array 18 is formed by 
detector elements 20 which together sense the project- 
ed x-rays that pass through a medical patient 22. Each 
detector element 20 produces an electrical signal that 
represents the intensity of an impinging x-ray beam and 
hence the attenuation of the beam as it passes through 
patient 22. During a scan to acquire x-ray projection da- 
ta, gantry 12 and the components mounted thereon ro- 
tate about a center of rotation 24. 
[0019] Rotation of gantry 12 and the operation of x- 
ray source 1 4 are governed by a control mechanism 26 
of CT system 10. Control mechanism 26 includes an x- 
ray controller 28 that provides power and timing signals 
to x-ray source 1 4 and a gantry motor controller 30 that 
controls the rotational speed and position of gantry 12. 
A data acquisition system (DAS) 32 in control mecha- 
nism 26 samples analog data from detector elements 
20 and converts the data to digital signals for subse- 
quent processing. An image reconstructor 34 receives 
sampled and digitized x-ray data from DAS 32 and per- 
forms high speed image reconstruction. The recon- 
structed image is applied as an input to a computer 36 
which stores the image in a mass storage device 38. 
[0020] Computer 36 also receives commands and 
scanning parameters from an operator via console 40 
that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the recon- 
structed image and other data from computer 36. The 
operator supplied commands and parameters are used 
by computer 36 to provide control signals and informa- 
tion to DAS 32, x-ray controller 28 and gantry motor con- 
troller 30. In addition, computer 36 operates a table mo- 
tor controller 44 which controls a motorized table 46 to 
position patient 22 in gantry 12. Particularly, table 46 
moves portions of patient 22 through gantry opening 48. 
[0021] The following discussion which describes cor- 
recting for z-axis non-uniformity sometimes refers spe- 
cifically to an axial scan. The correction algorithm, how- 
ever, is not limited to practice in connection with only 
axial scans, and may be used with other scans, such as 



helical scans. It should be further understood that the 
algorithm described below may be implemented in com- 
puter 36 and would process, for example, reconstructed 
image data. Alternatively, the algorithm could be imple- 

s mented in image reconstructor 34 and supply corrected 
image data to computer 36. Other alternative implemen- 
tations are, of course, possible. 
[0022] As shown in Figure 3, detector array 18 in- 
cludes a plurality of detector modules 20 secured to an 

10 arc shaped detector housing 50. Each detector module 
20 includes a multidimensional photodiode array (not 
shown) and a multidimensional scintillator array (not 
shown) positioned in front of and adjacent to the photo- 
diode array. The photodiode array includes a plurality of 

15 photodiodes (not shown) which are optically coupled to 
the scintillator array. Each photodiode, identified as a 
cell, generates an electrical output signal representative 
of the light output by each scintillator of the scintillator 
array. In one embodiment, the cell output signals are 

20 supplied through a field effect transistor (FET) array (not 
shown) to DAS 32. By controlling the FET array, specific 
photodiode outputs are electrically supplied to DAS 32. 
In one embodiment where detector array 1 8 includes fif- 
ty-seven detector modules 20, with each detector mod- 

25 ule 20 including a photodiode array having 16 columns 
per channel, the FET array is configured so that all 16 
output signals of the photodiode array, identified as a 
channel, are electrically supplied to DAS 32 to generate 
16 simultaneous slices of data. Of course, many other 

30 combinations are possible. 

[0023] For example, the FET array may be configured 
to provide other numbers of slices, including one, two, 
and four slice modes. Depending upon the specific con- 
figuration of the FET array various combinations of pho- 

35 todiode outputs may be enabled, disabled and com- 
bined so that the thickness of each slice may be 1 , 2, 3, 
or 4 rows. Additional examples include, a single slice 
mode including one slice with slices ranging from 1 row 
to 16 rows thick; and a two-slice mode including two slic- 

40 es with slices ranging from 1 row to 8 rows thick. Addi- 
tional modes beyond those described are possible 
where the total number of photodiode array element 
rows, or pixels per channel, is equal to the number of 
slices or FET array outputs times the number of rows 

45 per slice. For example, in a 4 slice mode of operation 
using 4 rows per slice, photodiode and scintillator arrays 
include at least 1 6 rows of cells and FET array includes 
at least four outputs. In one embodiment, for example, 
each row is 1 .25 mm wide. 

50 [0024] A calibration algorithm may be used to correct 
projection data for z-axis, or z-direction non-uniformity. 
Specifically, the calibration algorithm may be used to 
correct for z-axis non-uniformity of the x-ray beam and 
detector array 18. More specifically, utilizing the output 

55 signals of the photodiode array, separate z-axis non-uni- 
formity profiles are generated for the x-ray beam and 
detector array 18. 

[0025] Referring to Figure 4 and in one embodiment, 
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the calibration algorithm generates the z-axis x-ray 
beam non-uniformity profile by obtaining, or sampling 
the output signals from at least one channel in the z- 
direction of detector array 18. Specifically, prior to com- 
bining the output signals with the FET array as deter- 
mined by the configuration of detector array 1 8, the out- 
put signal from each detector cell of at least one chan- 
nel, referred to as a fully-sampled channel, of array 18 
is supplied to DAS 32. Utilizing the output signals from 
the fully sampled channel, a z-axis x-ray beam profile is 
determined for the channel. For example, where detec- 
tor array 18 is configured to include four slices having a 
slice thickness of 5mm (4 - 1.25mm rows), the output 
signals are supplied to DAS 32 from each of the 1 6 cells 
in at least one channel for determining the z-axis x-ray 
beam profile. After the output signals are sampled, the 
output signals are combined to generate the four slice 
5mm images. 

[0026] Using the output signals from a plurality of fully 
sampled channels, a z-axis profile may be generated for 
the entire x-ray beam. Specifically by fully sampling a 
plurality of channel output signals a profile of the incom- 
ing x-ray beam flux may be estimated using various 
curve fitting and interpolation methods. More specifical- 
ly because the incoming x-ray beam profile in the z-axis 
direction is a slowly varying function of the channel 
number the z-axis profile variation can be assumed to 
be almost constant for a small group of channels adja- 
cent to each fully sampled channel. Particularly and in 
one embodiment, the output signals from the fully sam- 
pled channels are interpolated to generate the x-ray 
beam z-axis profile for all non-fully sampled channels. 
Using the resulting x-ray beam z-axis profile, the output 
signal from each detector cell may be corrected for x- 
ray beam z-axis non-uniformity. In other embodiments, 
the x-ray beam z-axis profiles may be generated using 
curve fitting methods as well as assigning the fully sam- 
pled channel z-axis profile to non-fully sampled adjacent 
channels. 

[0027] In an alternative embodiment, less than all of 
the detector cells in a selected channel are sampled. 
Specifically, after selecting the channels to be sampled, 
at least one cell of each slice is sampled to generate the 
x-ray beam z-axis profile. More specifically and in one 
embodiment shown in Figure 5, where detector array 18 
is configured in a four slice 5mm thickness mode, prior 
to combining the output signals with the FET array, one 
output signal from each slice is supplied to DAS 32 for 
generating the x-ray beam z-axis profile. Particularly 
prior to combining the output signals of respective cells 
1 through 4, cells 5 through 8, cells 9-12, and cells 1 3 
- 16 to generate data for respective slices 1, 2, 3, and 
4, one output signal is sampled from respective cells 1 
through 4, cells 5 through 8, cells 9-12, and cells 13 - 
16. For example, output signals may be sampled from 
respective cells 2, 6, 11 and 15 to generate an x-ray 
beam z-axis profile for the sampled channel. For exam- 
ple, using a known curve fitting or technique and the out- 



put signals from cells 2, 6, 11 and 15, may result in an 
x-ray beam profile as shown in Figure 6. Thereafter, by 
sampling multiple channels and using curve fitting or in- 
terpolation methods as described above, the sampled 

5 channel z-axis profiles may be used to generate z-axis 
profiles for the non-sampled channels. 
[0028] After determining the x-ray beam z-axis profile, 
the output signals for the configured slice and slice thick- 
ness mode are calibrated or adjusted based on the z- 

10 axis profile. Specifically, the measured output signals 
from the non-fully sampled channels are calibrated to 
adjust the signals for the z-axis profile. In one embodi- 
ment and as shown in Figure 7, detector array 18 in- 
cludes sixteen rows of cells and is configured in a four 

15 slice 5mm thickness mode, where four 1 .25mm cells are 
combined to form one 5mm cell. More specifically and 
referring to slice 1, the gains for the fully sampled cells 
are 9o> Gh> 92> 93 and the variations of x-ray flux in the 
z-direction, or x-ray beam z-axis profile, are oc 0 /, ot 1 l, oc 2 l, 

20 <x 3 l. Corresponding parameters for non-fully sampled 
cells are g' 0 , g'-|, g'2, 83, o! 0 \\ a\\\ a' 2 l\ and a' 3 r. The 
parameters g' 0 , g'-|, g^, g'3 may be generated using 
known gain calibrations for the non-fully sampled chan- 
nels. As described above, oc' 0 , a'-,, a' 2 , and a' 3 may be 

25 generated using interpolation, curve fitting or may be as- 
signed the values of adjacent fully sampled cells. The 
measured output signal from a non-fully sampled chan- 
nel y 1 is: 

30 

M 
/ = 0 



l'= average x-ray intensity of a non-fully sampled 
cell, 

40 g'j = gain of non-fully sampled channel cell i, 
oc'j = x-ray beam z-profile variation on cell i. 

[0029] Specifically, and in one embodiment, detector 
array 18 is configured in a four slice mode, where four 
45 1.25mm cells are combined to form a 5mm slice, the 
measured output signal from a non-fully sampled chan- 
nel, y\ is: 



so y'= So^o'' + 0i a V + ffW + 03 a y' 

[0030] The true x-ray flux to the channel, x', is: 

55 



30 



35 

where: 
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M 




1=0 

[0031] In the four 5mm slice mode described above, 
the true x-ray flux to a non-fully sampled channel, x\ is: 

x , = a l 0 + a\ + a l 2 + a l 3 

Thetrue x-ray flux and measured output signals are sim- 
ilarly determined for other slices, for example, slices 2, 
3 and 4. 

[0032] As a result of the slowly varying z-profile vari- 
ation, the a parameters can be estimated on multiple 
views or an average of multiple views to reduce the in- 
fluence of statistical fluctuation. In addition, the use of 
multiple views reduces the computational burden of the 
correction algorithm. 

[0033] Various methods may be used to determine 
which cells may be sampled to determine the z-axis pro- 
files for the x-ray beam and detector array 18. In one 
embodiment where multiple channels are fully sampled, 
the z-axis uniformity of the channels may be utilized to 
determine the spacing of the fully sampled channels. 
Specifically, because the center channels are more sen- 
sitive to z-axis non-uniformity, the spacing of the chan- 
nels to be fully sampled are closer together near the iso- 
center channel and farther apart away from the isocent- 
er channel. In another embodiment, information from 
adjacent channels may be used to generate the z-axis 
profile so that the error introduced from any one cell is 
minimized. 

[0034] In an additional embodiment, only those detec- 
tor cells known to have uniform gain measurements are 
manually selected for generation of the z-axis profile. 
Specifically, because the output of these uniform gain 
channels have no error, the channel is not required to 
be fully sampled. Utilizing these cells, a 0 , a 1 , a 2 , and oc 3 
may be estimated based on a polynomial fit of the meas- 
ured signal from the uniform gain cells. After estimating 
the measured signal, interpolation, or one of the other 
methods described above, may be utilized to estimate 
a' 0 , a'-, , a' 2 , a' 3 for the non-uniform cells. 
[0035] In another additional embodiment, an estimate 
of the x-ray beam incoming flux profile is generated and 
used to calibrate the detector cells output signals. The 
calibrated measurements of the detector cells are then 
used to generate a revised x-ray beam z-axis profile. 
[0036] The above described system includes a cali- 
bration algorithm for correction of error due to z-axis 
non-uniformity of the x-ray beam and the detector array. 



In addition, the system enables correction of varying de- 
tector cell gain profiles. As a result, artifacts are reduced 
in a displayed image. 

[0037] From the preceding description of several em- 
5 bodiments of the present invention, it is evident that the 
objects of the invention are attained. Although the in- 
vention has been described and illustrated in detail, it is 
to be clearly understood that the same is intended by 
way of illustration and example only and is not to be tak- 
en by way of limitation. For example, the CT system de- 
scribed herein is a "third generation" system in which 
both the x-ray source and detector rotate with the gantry. 
The present invention, however, may be used with many 
other CT systems including "fourth generation" systems 
wherein the detector is a full-ring stationary detector and 
only the x-ray source rotates with the gantry. The 
present invention could also be utilized in connection 
with stop-and-shoot as well as helical scanning type CT 
systems. In addition, although the correction algorithm, 
in one form, has been described being performed on da- 
ta directly from the output signals, the present invention 
could be implemented at various points in the data cor- 
rection/processing process. 

Claims 

1. A system for calibrating projection data for generat- 
ing a tomographic image of an object, said system 
comprising an x-ray source and a detector array 
comprising a plurality of cells arranged in rows and 
columns and displaced along a z-axis, said system 
configured to: 

determine at least one z-axis x-ray beam pro- 
file; and 

determine a true x-ray flux utilizing detector 
gain profiles and said determined x-ray beam 
profile. 

2. A system in accordance with Claim 1 wherein said 
detector cell columns are arranged in channels and 
wherein to determine said x-ray beam profile, said 
system is configured to obtain signals from a select 
number of detector cells in at least one channel of 
said detector array. 

3. A system in accordance with Claim 2 wherein each 
channel comprises at least one slice and wherein 
to obtain signals from a plurality of detector cells in 
at least one channel, said system is configured to 
obtain a signal from at least one cell in each slice. 

4. A system in accordance with Claim 3 wherein to ob- 
tain a signal from at least one detector cell in each 
slice, said system is configured to obtain a signal 
from each of said detector cells in said channel. 
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5. A system in accordance with Claim 2 wherein to de- 
termine detector gain profiles, said system is con- 
figured to determine a gain profile for each of said 
selected detector cells. 

6. A system in accordance with Claim 1 further config- 
ured to adjust the projection data using said deter- 
mined detector gain profiles and said x-ray flux. 

7. A system in accordance with Claim 1 wherein to ad- 
just the projection data, said system is configured 
to generate a gain profile for each nonselected de- 
tector cell. 

8. A system in accordance with Claim 7 wherein to 
generate a gain profile for each non-selected detec- 
tor cell, said system is configured to determine a 
measured output signal, wherein said measured 
output signal, y' ; is: 

y=f g'ia'j' 

where: 



1 2. A system in accordance with Claim 1 0 wherein said 
true x-ray flux for said channel is: 

5 x >= ° ] £ 3 , 

S^O + 0'l a 'l + 9 l 2 a '2 + 9 l 3 al 3 

13. A system in accordance with Claim 12 wherein to 
generate said xray flux profile for said non-selected 

10 detector cells, said system is configured to create a 
polynomial fit of said selected detector cell output 
signals. 

1 4. A system in accordance with Claim 2 wherein to ob- 
15 tain signals from a selected number of detector 

cells, said system is further configured to: generate 
view data utilizing said detector cell signals; and 
generate said true x-ray flux profile utilizing said 
view data. 

20 

1 5. A system in accordance with Claim 2 wherein to ob- 
tain signals from a selected number of detector 
cells, said system is further configured to: generate 
data for multiple views utilizing said detector cell 

25 signals; and generate said true x-ray flux profile uti- 
lizing said multiple view data. 



g'j = gain of cell i in detector channel; 

a'j = variation of flux in z direction on cell i; 

I 1 = average x-ray intensity of a non-fully sam- 30 

pled cell; and 

M = number of cells that form a slice. 



1 6. A system in accordance with Claim 2 wherein to ob- 
tain signals from a selected number of detector 
cells, said system is further configured to select at 
least one detector cell having an uniform gain 
measurement. 



9. A system in accordance with Claim 8 wherein to 
generate said true x-ray flux for each non-selected 35 
detector cell, said system is configured to determine 

a true x-ray flux for said channel. 

10. A system in accordance with Claim 9 wherein said 
true x-ray flux for said channel is: 40 



1 7. A system in accordance with Claim 2 wherein to de- 
termine an x ray beam z-axis profile, said system is 
configured to: 

estimate an x-ray flux profile; 
adjust said output signals of each channel; and 
generate a revised x-ray beam z-axis profile uti- 
lizing said adjusted output signals. 



M 
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1 1 . A system in accordance with Claim 1 0 wherein each 
slice channel comprises four cells and wherein said 
measured output signal for each slice of each non- 
selected detector cell is: 



18. A method for calibrating projection data in a multi- 
slice computed tomography system, the system in- 
45 eluding an x-ray source and a detector array com- 
prising a plurality of cells arranged in rows and col- 
umns and displaced along a z-axis, said method 
comprising the steps of: 

50 determining an x-ray beam z-axis profile; and 

determining a true x-ray flux utilizing detector 
gain profiles and the determined x-ray beam 
profile. 

55 19. A method in accordance with Claim 1 wherein the 
detector array column cells are arranged in chan- 
nels and wherein determining thex-ray beam profile 
comprises the step of obtaining signals from a se- 
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lect number of detector cells in at least one channel 
ot the detector array. 

20. A method in accordance with Claim 1 further com- 
prising the step of adjusting the projection data us- s 
ing the determined detector gain profiles and xray 
flux profiles. 
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